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Abstract. We introduce a nonlinear structured population model with diffusion in 
the state space. Individuals are structured with respect to a continuous variable 
which represents a pathogen load. The class of uninfected individuals constitutes a 
special compartment that carries mass, hence the model is equipped with generalized 
Wentzell (or dynamic) boundary conditions. Our model is intended to describe the 
spread of infection of a vertically transmitted disease, for example Wolbachia in a 
mosquito population. Therefore the (infinite dimensional) nonlinearity arises in the 
recruitment term. First we establish global existence of solutions and the Principle 
of Linearised Stability for our model. Then, in our main result, wc formulate simple 
conditions, which guarantee the existence of non-trivial steady states of the model. 
Our method utilizes an operator theoretic framework combined with a fixed point 
approach. Finally, in the last section we establish a sufficient condition for the local 
asymptotic stability of the positive steady state. 



1. Introduction of the model 

Structured population dynamics is an exciting field of research in mathematical bio- 
logy, see e.g. [T2j [29j [32l [38] for a collection of classical models and results in the area. 
Structured population models often assume spatial homogeneity of the population in a 
given habitat and only focus on the dynamics of the population arising from differences 
between individuals with respect to some physiological structure. In a general model, 
reproduction, death and growth characterize individual behaviour, which may be affected 
by competition (scramble or contest), for example for available resources or for a mating 
partner. One usually incorporates certain types of nonlinearities via density dependence 
in the vital rates to account for these biological phenomena. 

Traditionally, structured population models have been formulated as first order non- 
linear partial differential equations of hyperbolic type and therefore the analysis of these 
models is often challenging from the mathematical point of view. In the very recent 
paper |27j Hadeler introduced size-structured population models with diffusion in the 
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size-space. The biological motivation is that diffusion allows for "stochastic noise" to be 
incorporated in the model equations in a deterministic fashion. In |27j Hadeler mainly ad- 
dressed the question that what type of boundary conditions are necessary to be imposed 
for a biologically plausible and mathematically sound model. In this context particular 
cases of a Robin boundary condition were considered. We note that other researchers 
also investigated recently the effects of introducing stochasticity (for example variation 
in individual growth rates) via diffusion in structured population models, see 

In this paper we introduce and analyze a structured population model, with so called 
distributed recruitment term and generalized Wentzell boundary conditions. Our model 
intends to describe the dynamics of a population, which is infected with a certain type of 
bacteria, see later for more details. At the same time our model is general enough to make 
the forthcoming analysis to be applicable for similar classes of models. In particular, we 
consider the following model: 

ut{x, t) + (7(a;)w(x, t))^ = {d{x)u^{x, t))^ - n{x)u{x, t) 



+ /3o (Uit)) /3i(x, y)P2 (^^^77(^) j ^) ^ e (0, m), 

(1.1) 

mx)u^{x, t))X^o - (rf(0) + 7(0))i^.(0, t) 

+ (^7(0) - y (0) - mt)) b2 (^^^4|r^) ) ^ °' (^-^^ 

d{m)ux{m, t) — j{m)u{m, t) ~ 0, (1-3) 

with a suitable initial condition. The function u = u{x, t) denotes the density of indi- 
viduals of pathogen load x at time t. This means that the total number of individuals 
with pathogen (or bacterium) load between Xi and X2 is given by J^^ u{y, t) dy. It is 
then clear that the uninfected individuals constitute a separate special class. Therefore, 
in contrast to any other state, the state a; = carries mass. Hence we specify the value 
of the function u at a; = 0, which gives the total number of individuals in the uninfected 
class. Then the total population size at any time t is given by 

U{t) = |m(0,OI+ / \u{r,t)\ dr = Hull;,. 

In mathematical terms, we distinguish between equivalence classes of the state space 
i^, which are represented by functions, which have different values at a; = 0. Or, in 
other words, to each equivalence class of the Lebesgue space we assign a real number, 
which is its "value" at a; = 0. This is expressed in the Wentzell- type boundary condition 
()1.2p . We also assume a maximal infection load m. The function 7 represents the 
reproduction rate of the bacterium in the host. Since the maximum infection load is m 
we have the boundary condition (|1.3p . i.e. the infection load cannot increase above m 
in any individual in the population. It is shown by a straightforward calculation (see 
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e.g. [13 for a similar calculation) that the boundary conditions (|1.2p - (jl.3|) guarantee 
conservation of the total population size in the absence of mortality and recruitment. As 
usual, ^ denotes the mortality rate (both natural mortality and extra mortality due to 
the infection) of individuals, while d stands for diffusion. Now wc impose the following 
general assumptions on the model ingredients: 

/3o,6oeC([0,oo)),/32,fe2eC([-l,l]),A*eC([0,TO]), /3i eC([0,m] x [0,m]), 
/3o,/3i,/32,6o,&2,A* > 0, 7,deCH[0,m]), d>0. 

In addition we will make some technical assumptions on some of the model ingredients 
later on. Our goal is to describe a model of a population, which is infected with a disease 
that is only transmitted vertically; hence the only nonlincarity arises in the term specify- 
ing the recruitment of individuals. To be more precise, we are interested in the dynamics 
of Wolbachia infections in mosquito populations. Wolbachia is a maternally transmitted 
endosymbiont bacterium, which infects around 70% of the arthropod species. Recent 
laboratory results suggest, see [21], that a stable introduction of a particular strain of 
Wolbachia into A. Aegypti halved the average life span of an adult mosquito. This repre- 
sents a potential tool to eradicate mosquito born diseases (like malaria, dengue or West 
Nile virus), since the life span of the mosquito is reduced to 2-3 weeks which more or 
less equals the intrinsic excubation period of the disease, i.e. the time needed for a 
mosquito to become infectious after feeding on an infected human. According to [5T] . 
high maternal inheritance, strong cytoplasmic incompatibility and low costs to reproduc- 
tive output arc the key factors for prevalence of the Wolbachia infection. Cytoplasmic 
incompatibility (CI for short) means that from a mating of an uninfected female and 
infected male there is no (complete CI) or only a very few (partial CI) number of viable 
offspring. We refer the interested reader to [HI dni [M] for more biological background 
and mathematical modelling of populations infected with a cytoplasmic incompatibility 
inducing Wolbachia. The possible underlying mechanisms for complete or partial expres- 
sion of CI are still a matter of debate for different species. In our model we take the 
following view of partial CI, see [SI Hi] for further reference. We make the assumption 
that a female who has bacterium load y can successfully mate only with males of bac- 
terium load less than y. Therefore the function /32, which determines the fertility rate of 
an individual with infection load y, depends on the proportion of the individuals, which 
have infection load higher than y. In particular /32 is a monotone decreasing function of 
its argument. In other words, we assume that the sex ratio is 1 : 1 and the probability for 
a female mosquito of producing a viable offspring is a monotone decreasing function of 
the proportion of the male population size that have higher infection load. The function 
/3q represents competition effects due for example to limitations in available resources. 
Similarly, 62 measures the reproductive success of uninfected females and 60 represents 
the effects of competition on the reproduction of uninfected individuals. The maternal 
transmission rate is determined by the function /3i, i.e. infected mothers with bacterium 
load y give birth to offspring with bacterium load x at a rate of /3i{x, y). 
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The boundary condition (|1.2|) can be rewritten, and more easily understood from the 
point of view of its biological meaning, in the dynamic form (by formally replacing the 
diffusion operator from equation (jl.ip on the boundary) as: 



It is not difficult to verify that in the absence of mortality and recruitment, i.e. taking 
/i = and bo = f3o = 0, the total population size U{t) is conserved. Notice that in 
this case, ut{0,t) equals — 7(0)u(0,i) + d{0)urc{0,t), i.e., the flux through a; = due 
to reproduction of the bacteria (from the biological point of view it would be natural 
to assume 7(0) = 0, but we do not impose such a restriction here) and to diffusion 
(noise) in the second term (this could take into account a certain -very weak- horizontal 
transmission) . 

Our main goal in this paper is to establish sufficient conditions for the existence 
of positive steady state solutions of model p.ip - (|1.3p . We shall refer the interested 
reader to [51[51 ll5[|17j where different size-structured models with distributed recruitment 
processes were investigated. The boundary condition (jl.2p is the so called generalized 
Wentzell-Robin or dynamic boundary condition. These "unusual" boundary conditions 
were investigated recently for models describing physical processes such as diffusion and 
wave propagation, see e.g. [HlUnillS]- Briefly, they are used to model processes where 
particles reaching the boundary of a domain can be either reflected from the boundary or 
they can be absorbed and then released after some time. Recently in |17| we introduced, 
as far as we know for the first time, Wentzell-type boundary conditions in the context of 
physiologically structured populations with a distributed recruitment process and with 
diffusion in the size space. The introduction of diffusion in the size space is very natural 
from the application point of view, since in the real world individuals who start their life 
in the same cohort do not finish their life so, due for example to stochastic variations in 
individual growth rates. We refer the reader to [37] (see also [5]) where different types 
of population models were introduced with diffusion in the size-space. 

The idea of considering population models where the structuring variable represents 
a pathogen load is clearly not new. In fact in [57] Waldstadter et al. introduced a 
similar model, where they have derived a dynamic boundary condition for the special 
class of uninfected individuals. In their model this class however does not correspond to 
the value u(0,t), in fact the uninfected population size is represented by a new variable 
U{t) and a transition condition between the compartment U and u(0,t) is given. This 
implies two things: firstly u{0,t) ^ U{t), secondly individuals who are sitting in the 
special compartment U are not subject to diffusion, i.e. first they need to enter the state 
u(0,t). This in some sense may seem counterintuitive, as diffusion is incorporated to 




(1.4) 
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model stochastic noise, which could result in low probability random infections. Since we 
model an (almost) completely vertically transmitted disease, we do not take into account 
infection in a usual way. By usual way we mean by means of a nonlinear transition 
rate of the type SI, which in our model would correspond to u{0, t) J™ u{x, t) Ax. But, 
due to the nature of the Wentzell boundary condition our model allows for transition 
between the uninfected and infected compartments as a diffusive flux through a; = in 
both directions (one of them corresponding to random infections and the other one to 
spontaneous, i.e. random healings), depending on the sign of u^, as commented above. 
We note that the mathematical analysis presented in [37| is carried out in the Hilbert 
space L^. 

We shall mention here that the first papers introducing boundary conditions that 
involve second order derivatives for parabolic or elliptic differential operators go back to 
the 1950s, see the papers by Feller [351 123] and Wentzell (also transliterated as "Ventcel' 
", Aleksandr D. Ventcelb) [3SJ[3S]. These first studies were purely motivated from the 
mathematical point of view. The original question, as far as we know, was to identify 
the maximal set of possible boundary conditions that give rise for a parabolic differential 
operator to generate a contraction semigroup on an appropriate state space. 

In this paper first we establish global existence and positivity of solutions of model 
(|l.ip - (|1.3p . Our existence proof relies on the existence of the semigroup governing the 
linear part of the model, which was established in |17j (with two point Wentzell boundary 
conditions) following similar arguments developed in [121 HOji and on the analyticity 
of this linear semigroup, which was very recently established in |21j (and in fact for 
higher dimensional domains, as well). Then in Section [3] we investigate the existence 
of non-trivial steady states to our model. To treat the steady state problem, we devise 
Schauder's fixed point theorem combined with an operator theoretic approach. A similar 
operator theoretic framework was previously utilised for simpler problems (in particular 
with finite dimensional nonlinearities and classical boundary conditions), see e.g. [5lll5j. 
The key idea to treat the steady state problem is to define a linear operator for a fixed 
environment (nonlinearity) and to study spectral properties of that operator. Finally in 
the last section we establish some sufficient stability conditions for the positive steady 
state. 



In this section we establish global existence (and positivity) of solutions of the nonli- 
near problem (|l.ip - (jl.3D . Throughout the section we employ some standard results from 
[28] . We introduce the state space X = L^(0, m) © R with norm 




2. Existence and positivity of solutions 




|u|| = ||(^,uo)| 



kill 4- |wo| 
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which is a Banach lattice. Next we write our problem (jl.l[) - ()1.3|) in the form of an 
abstract Cauchy problem: 

^-Au = F(u), t>0, u(0) = uo, (2.5) 
where A is the linear operator defined by 

" {-fi{0)u{0) + d{0)u'{0) - 7(0)u(0)) ■ ^^-^^ 
In (|2.6p the operator A is defined as 

The domain of the linear part A is given by 

D{A) = {u e C2[0,m] : Au G L^(0,m), d{m)u'{m) - 7(m)u(m) = 
idis)u'{s)y\^^^ - (rf(0) + 7(0))ii'(0) + (7(0) - I'muiO) = 0} . 

The nonlinear but bounded function F is defined as 

^w-C'"'".! ft "'^^'°'}. p**) 

with D(i^) = X, where 

/3o(||u||)/;'/3i(,y)/32(^^^^^^)«(y)dy 
^/?o(||u||)/™/3i(0,y)/32(^^^)«(2/)d2/ + 6o(||u||)^ 

Note that is only defined on the open set X \ {0}. 



Theorem 2.1. Assume that /3o,/32,^o '"^'^ ^2 o-fe locally Lipschitzian. Then a unique 
solution of problem (|2.5p exists for all positive times for any uq G X. 

Proof. We use [551 Corollary 3.3.5] to establish global existence of solutions to the 
semilinear problem (j2.5p . To this end, first we show that F is locally Lipschitz continuous, 
i.e. for every u € X there exists a neighbourhood U C of u such that for every 
Ui,U2 G [/ we have ||-F(ui) - F(u2)||;^. < L ||ui - U2||;^. = L (||ui - U2II1 + - W2ol) 
for some constant L > 0. 

If either Ui = or U2 = then since F is bounded (i.e. all the ingredients /3q, Pi, f32, &o 
and 62 are continuous, hence bounded) and F{0) = we have the Lipschitz property 
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||i^(u)|| < i||u|| for some L > 0. For Ui 7^ 0,U2 7^ we have: 
||F(ui)-F(u2)|U < 



|/3o(||ui||)-/3o(||u2| 



+ /?0(||U2||) 



ui{y) - P: 



ui{y) dy 



|U2| 



+ |/3o(||ui||)-/3o(||u2||)| 



+ /3o(||u2| 



/3i(0,2/) /32 



/5i(0,y)/32 ^- 
/ /""iW dr 



/J"Hi(r) dr \ 
lluill ) 



ui{y) dy 



"2 



+ |5o(||ui||)-6o(||u2||)|62 

/J"m(r) dr 



&0(||U21|) 62 



/J"Mi(r) dr ' 
lluill 

Ml,, - 62 



/"M2(r) dr ' 
IIU2II 



U2o 



(2.9) 

"2(2/)^ dy 
(2.10) 

(2.11) 

W2(y)^ dy 
(2.12) 

(2.13) 
(2.14) 



The term in (|2.9p can be bounded above by i^(,||ui — U2||Af^jAf/32m||Mi|Ii, where /?o < 



M^Q,/3i < M/3j^,/32 < Af^2 and is the Lipsehitz constant of /3o in 
obtain the estimate for (|2.10p we note that 



0, sup ||m 
ueu 



To 



^2 



7™wi(r)dr' 



/32 



|Ul|| ^ 

7"^i(0 dr ' 
lluill 



wi(y) - ^2 



7"M2(r) dr' 



/S2 



^ IIU2II 
7" "2(r) dr ' 
IIU2II 



M2(y) 



ui{y) + (32 



7J""2(r) dr' 



|U2| 



(ui(y) - ^2(2/)) 



<Lf, 



/"Mi(r)dr /™U2(r)dr /"u2(r)dr /™M2(r)dr 



Ul 



Ul 



Ul 



U2 



\ui{y)\ +MfiJui{y) - W2(y)| 



+ ^/32 



Mi(r)dr— I U2{r)dr 
y 

/J"u2(r) dr 



U2 - Ul 



|U2| 



\ui{y)\ \ 

1 lluill ; 

1 \My)\ 

" lluill 



+ Mp^luiiy) - M2(y)| 



(2.15) 
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where is the Lipschitz constant of P2 in [—1,1]. Similarly to (j2.15p . we obtain the 
appropriate estimates for the terms (|2.1ip - (|2.14p . Hence F is locally Lipschitzian. Since 
|F(u) — -F(0)|| — ||F(u)|| < i||u||, F is also sublinear and since A is sectorial (see pT] ) 
the statement of Theorem 12. II follows from [551 Corollary 3.3.5]. □ 
Note that F : The solution of the Cauchy problem (j2.5p can be written 

as 

u{t)=T{t)uo+ f T{t-s)F{u{s))ds, iG(0,to), (2.16) 
Jo 

where T is the linear semigroup generated by the closure of the sectorial operator A. 
We refer the reader to |17| were it is showed that a similar linear problem (with two 
point Wentzell boundary condition) is governed by a quasicontractive positive semigroup. 
Furthermore we refer to |21j where the sectoriality of the operator A was shown in 
L^, 1 < p < 00, in general. Since is a positive operator and the semigroup T{t) is 
positive the variation formula (|2.16p immediately shows positivity of solutions. 

The Principle of Linearised Stability can be established for model (|l.ip - (|1.3p by di- 
rectly applying again results from [55] ■ In particular, if is an equilibrium point then 
let F(^^ denote the linearisation of F at the equilibrium u*, which is a bounded linear 
operator defined on X (if it exists). See later in Section 4 for more details. Then [28l 
Theorem 5.1.1] implies that if A is sectorial and F is locally Lipschitz continuous then 
the equilibrium u, is asymptotically stable if <t (A + F^J C {A G C | i?e(A) < a < 0}. 
On the other hand, [551 Theorem 5.1.3] implies (under the same conditions on A and F) 
that the equilibrium u, is unstable if cr (A + F(^^) n {A G C | i?e(A) > 0} is not empty. 



3. Existence of non-trivial steady states 

It is obvious that model (|l.ip - (|1.3p admits the trivial steady state. It is also clear that 
even the time independent version of equations (|l.ip - (|1.3p cannot be solved explicitly. 
Therefore, to establish conditions which guarantee the existence of a positive steady 
state we utilise a combination of an operator theoretic framework (see e.g. [3 [151) 
a fixed point approach. For basic definitions and results from linear semigroup theory 
used throughout this section we refer the reader to [H [TUl [T3] . 

For a fixed non vanishing v G L^(0,m) © M let us define the linear operator 4'v 
(parametrised by v) by 

«'vU = Au + F°u, (3.17) 
where the operator A is defined in (j2.6p . and 



F2n^ 



,/3o(||v||)/™/3i(0,,)/3,(^^::^).(y)d,^ 
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with domain 

D{-^^) ^ {u e C'^{[0,m]) : Au <E L^{0,m), d(m)u'(m) - 7(77i)u(m) = 

{dis)u'{s))X^^ (d(0) + 7(0))«'(0) + (7(0) - 7'(0))^.(0) = 0} . (3.18) 

It is clear that (u*,f*(0)) is a nontrivial steady state of model (|l.ip - (jl.3p if (and only if) 
it is a positive eigenvector belonging to the kernel of the (closure of the) linear operator 

Next we establish some (necessary) regularity properties of the semigroup generated 
by the closure of ^'v Several characterisations of irreducibility of a positive semigroup 
exist in the literature, see for example in O [101 US] • Here we follow [13] . A positive 
semigroup T on the Banach lattice X is said to be irreducible if the resolvent of its 
generator G is strictly positive, i.e. VO ^ / G X+ we have {R{X,G)f){x) > for almost 
all X and some A > s{G). 

Lemma 3.2. For every fixed v G Af-f the closure of the linear operator generates an 
irreducible semigroup on X . 

Proof. We introduce the mortality operator M as follows: 

^"=(-M0M0)) (3-19) 

We then consider the resolvent equation 

(AI - (A - A/)) u = h, (3.20) 

for h G X+, A > 0, i.e. 

—h{x)~{d{x)ux(x)~^(x)u{x)),^~\u{x), xg(0,to), (3-21) 
-ho = - (7(0) + A) w(0) + d(0)w'(0), (3.22) 

for an unknown u G D{A). We have shown in |17l Theorem 2.1] that the linear semigroup 
generated by A — M is positive (with two point Wentzell boundary conditions, but the 
proof can be trivially adapted to our case), i.e. the resolvent operator {XL — {A — M))~ 
is positive for A > large enough. Hence the solution u of equations p.2ip - p.22l) is 
non-negative, i.e. u > 0. It is only left to show that the solution u is in fact strictly 
positive. 

The one-dimensional minimum principle assures that u cannot attain its minimum 
value at an interior point of (0,m), see e.g. [33l Ch. 1 Theorem 3]. On the other 
hand, |33l Ch. 1 Theorem 4] guarantees that if u attains its non-positive minimum at 
then u'ifS) > holds, which contradicts h > 0, whereas the boundary condition at 
X = m implies u{m) > by the same theorem. Therefore the semigroup generated by 
the closure of A — M is irreducible. Since M is a bounded multiplication operator and 
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F° is bounded and positive for every v e X+ it follows that the closure of generates 
an irreducible semigroup for every v e see e.g. [51 C-III Proposition 3.3]. □ 

Remark 3.3 First of all we note that every steady state {v, v{0)) of model ()l.ip - (|1.3p 
shall have regularity W'^'^{0,m). This immediately implies that model (|l.ll) - p.3p does 
not admit a steady state which has uninfected individuals only, i.e. of the form (0, v(0)), 
with v{0) + 0. 

Secondly, as we noted before, a vector v ~ (u,u(0)) is a steady state of p.ip - p.3p 
if and only if it is a positive eigenvector belonging to the kernel of the closure of 
Therefore, irreducibility of the semigroup generated by the closure of shows that we 
also cannot have a steady state of the form (i;, 0), since every positive eigenvector of the 
closure of the generator is strictly positive if the semigroup is irreducible. 

Lemma 3.4. For every non vanishing v G the spectrum of 'fv can contain only 
isolated eigenvalues of finite algebraic multiplicity. 

Proof. Since M and F° are bounded it is enough to show that i?(A, A — M) is compact. 
This follows however by noting that the solution of the resolvent equation (|3.2ip - (|3.22p 
is in VF'^'^(0,to) ©R which is compactly embedded in X. The statement now follows on 
the grounds of [HI Proposition II.4.25] and [H Corollary IV.1.19]. □ 
As we noted before our goal is to show that there exists an element v G Xj^ such that 
the operator has eigenvalue zero. Then, Lemmas 13.21 and 13.41 guarantee the existence 
of a corresponding strictly positive (unique normalized) eigenvector. In case of a model 
with one dimensional nonlinearity, such as the one we treated in jl5j , the operator is in 
fact parametrized by a scalar quantity, namely, the total population size. In this case the 
positive steady state is obtained readily after an appropriate normalization of the positive 
eigenvector. In case of the model treated here the function /32 naturally depends on an 
infinite dimensional variable. Therefore, in general there is no guarantee that the positive 
eigenvector corresponding to the zero eigenvalue of the operator is in fact v (or a 
scalar multiple of it). For this reason, we need to construct an appropriate nonlinear map 
on a certain level set of the positive cone of the state space X and establish the existence 
of a fixed point of this map. We note that a different fixed point strategy was employed 
very recently in [18] for a class of models with infinite dimensional nonlinearitics (and 
zero flux boundary condition). That method, which uses fixed point results of nonlinear 
maps in conical shells of Banach spaces, does not apply to our model here since the 
construction of the nonlinear map requires the (implicit) solution of the steady state 
equation. We also note that a similar argument to the one used here, even though in the 
case of finite dimensional interaction variable is used in [7]. 

Theorem 3.5. Assume that Pq <^i^d, bo are strictly monotone decreasing functions of 
their argument and 

(i) lim I3q{x) = lim 5o(x) = and ^{x) > /iq > 0, Va; G (0,m); 
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(ii) there exists an r > such that for all G A'-f, ||v»|| < r we have that the 
spectral bound s > 0. 

Then model (|l.ip - (jl.3p has at least one non-trivial steady state. 

Proof. Since ^'v is a generator of a positive and irreducible semigroup with compact 
resolvent it follows that the spectrum cr(\I'v) is not empty, see e.g. [U C-III Theorem 3.7]. 
Moreover, the spectral bound s(^'v) is an isolated eigenvalue of algebraic multiplicity one 
with a corresponding strictly positive eigenvector and it is the unique eigenvalue with 
positive eigenvector, see e.g. [TUl Theorem 9.10]. We also note that the spectral bound 
and the corresponding positive eigenvector change continuously with respect to (the 
parameter) v, see e.g. [SHI Sect. 3 in Ch.4]. 
We introduce the level set 



It is shown that the closure of the linear operator A—M generates a (positive) contraction 
semigroup, hence we have s{A — RI) = ljo < 0. Intuitively it is clear that contractivity 
follows simply because in the absence of mortality and recruitment the total population 
size is preserved. Therefore, if 



then it is shown that s (^'v) < 0. Hence conditions (i) imply that there exists an i? > 
such that for Vv* G X+, ||v*|| > i? we have s (^'v.) < 0. From (ii) it then follows that 
5 C {x G \ r < ||x|| < R} . It also follows from conditions (i) and (ii) and from the 
continuity and strict monotonicity of the spectral bound (see below) that along every 
positive ray TZ = {av* | a G M+, G \ {0}} there exists a (unique!) v such that 
s('I'v) = 0. So the set S intersects every positive ray 7?. in a unique element. 

Next, we shall show that the spectral bound is in fact strictly monotone decreasing 
along every positive ray. To this end let < ai < a2 be real numbers and let v G X+. 
Consider the operators ^'a^v and ^'ajv- Both of them have compact resolvents by Lemma 
13.41 and generate positive and irreducible semigroups. We also note that ^'^^v ~ ^a2v is 
a positive operator, since 



^={xgA'+|s(*,,) = 0}. 




(/5o(ai||v||)-/3o(a2||v||))/™/3i(.,y)/32(. 
(/3o(ai||v||)-/3o(a2||v||))/7/3i(0,y)/32( 

( ' 

l^(6o(«l||v||)-&o(«2||v||))&2 (^S^) 



/J" "('') dr 



) 



u{y) dy 
u{y) dy 



) 



V 



l|v|| 





>o. 
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for every u G X^. We also have 

for A large enough. Hence [U Propositfon A. 2] implies that s(4'qiv) > s(4'q2v)- 

Next we construct the nonlinear map $ : S S, $(v) = {f^ o /a o /2 o /i)(v) = v', 
as illustrated briefly in the following diagram: 

$ : ^ *v Vv Vv ^-JC^- (3-23) 

eL(n(*v),A') eWl-\0,m)(BM eL^(0,m)©K 

For every element v £ 5 the map /i assigns the corresponding linear operator 
"parametrized" by v, which is defined via p.f 7p . This map /i is continuous and bounded. 
The map /2 assigns the strictly positive normalized eigenvector corresponding to the zero 
eigenvalue of the operator \l/v. This map /2 is clearly bounded and it is also continuous, 
in fact it is even analytic, see Lemma 1.3]. The map /s is the compact injection of 
W^'^{0, m) © R into L^{0, m) © R. Finally /4 is the projection along positive rays of the 
eigenvector Vv back into the set S, which is again continuous and bounded. 

Next we apply Schauder's fixed point theorem to an appropriately defined map. As 
we noted above the map $ is continuous and compact. It is only left to show that 
S is homeomorphic to a convex set. To this end, we define the map h : S B hy 
h{\i) = ji^, where B is the unit sphere intersected with the positive cone X^, i.e. 

S {u = {u,u^) e X+ I ||u|| = ||u||i + \u°\ = 1} . 

Then h is clearly continuous and one to one, since the spectral bound s (^') is strictly 
monotone decreasing along positive rays in X and every positive ray of X intersects B 
in a unique element. 

We shall show now that the function : B ^ S defined via /i^^(w) = aw, 

a G M+ such that s {^aw) ~ 0, is also continuous. Let Wn £ B he a. sequence such that 
Wn w G B and consider the sequence /i~^(w„) = q;„w„. It follows from condition (i) 
that there exists an i? > such that s(\['v) < for every ||v|| > R, hence the sequence 
a„ is bounded. Let be a convergent subsequence of and let a„j. — a. Since 
s ^^'q„^w„^^ = for every fc G N it follows from the continuity of the spectral bound 

that s (^Qw) = 0. Since there is exactly one element on the positive ray spanned by 
w at which the spectral bound vanishes we have d = a. If a„j is another convergent 
subsequence of a„ then the continuity of the spectral bound implies again that — > a. 
So a„ — a and finally /i~^(w„) = anW„ — 7> aw = /i~^(w). 

We note that the set B is convex since we are in an AL-space, i.e. we have ||f + g|| = 
||f|| + ||g|| for every f , g G X^^, sec e.g. [5]. Finally we apply Schauder's fixed point 
theorem (see e.g. [24]) to the continuous and compact map ^ : B ^ B defined by 

$(x) = /io$o/i-i(x). 
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to obtain a fixed point G B of the map which yields a fixed point = (x*) 
of the map $ in S. This x* is the positive steady state of model (|l.ip - (jl.3|) (see the 
characterization of a nontrivial steady state immediately after equation (j3.18p V □ 

Remark 3.6 Condition (i) is natural from the biological point of view. The first 
condition in (i) requires that the fertility rate of both infected and uninfected individuals 
tends to zero as the population size tends to infinity. This may be due for example to 
competition effects. In fact it turns out that as expected, cytoplasmic incompatibility 
itself does not have a negative feedback on population growth. The assumption of a 
strictly positive mortality function in (i) seems also realistic. 

Condition (ii) seems to be also a natural one, if one is to expect the existence of a 
positive steady state. In fact, if it is not satisfied then one can show that there exists 
a monotone decreasing sequence of positive real numbers r„ — > 0, such that for every 
n G N there exists a u" with ||u"|| < r„ and s (^E'uj) < holds. Then, the continuity 
of the spectral bound implies that s {^o) < and since the spectral bound is a strictly 
monotone decreasing function along positive rays, we have s (^'v) ^ 0, for every v e Af-f. 



4. Stability 

In the previous section we established conditions which guarantee the existence of a 
non-trivial steady state. The next natural step is to study the stability of the steady 
state (and also the stability of the trivial steady state). In this section we are going to 
establish a sufficient condition for the local asymptotic stability of the non-trivial steady 
state. First we note that the operator F is not (Frechet) differentiable at because it 
is homogeneous of degree 1 (it is an "angular operator"). Nevertheless it is Gateaux- 
diffcrentiable at 0. In particular its directional derivative at into the direction of any 
element v in the positive cone is: 

^( /3o(0)r/3i(.,y)/3.(^^^)%)d, 

°" [mo) /o" /3i(0, y)P2 {^^^) v{y) dy + 60(0)62 (^'^^) -0 

which is a nonlinear operator. 

If is a non-trivial equilibrium of (|l.ip - (jl.3p then we can formally linearise equation 
(jl.ip around . In particular if we denote by the linearisation (the Frechet deriva- 
tive) of the nonlinear operator F at then the linearised problem can be cast in the 
form of an abstract Cauchy problem 

^v=(yl + F^Jv, v(0)=vo, (4.24) 

where A is defined in (|2.6p and Z?(F^_) ~ X. We leave it for future work to address 
the linear problem (|4.24p (i.e. to address stability questions of equilibria), in general, 
for example using the Liapunov function techniques elaborated in [38] Sect. 4.2]. Here 
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however, we establish a straightforward and simple condition which guarantees that the 
non-trivial equilibrium is locally asymptotically stable. 

To this end wc consider on the state space X the operator F defined in (|2.8p with 
domain D{F) — . Then using approximations (in the spirit of [3S1 Dcf.2.4 in Sect. 
2.6]) such as 



/3o(||u,+v||) ^/3o(||u,||)+/3^(||u,||) / v{y)dy + v^ 



and similarly 

^ / /" ('^) +v{r)dr\ f /" (r) dr " 



/" Mr) dr\ f /" vir) dr p ^ ^J^^\iy)dy + vo^ 



it is shown that the linearisation of the operator F at the non-trivial equilibrium u.^ 
(according again to [38l Def. 2.4 in Sect. 2.6] is 

'rM-:y)P2{^^^i^) {M\K\\Hy) + mM\)Uo^viy)dy + vo)u4y)) dy' 

JJ"/3i(0,;/)/32( ^";^;;i'"' ) (/3o(||u,||Hy)+/?^(||u,||)(/;%(y)dy + «o)^*(y)) dy^ 
i(-,?/)/?o(l|u.||)/3^(^;^) {^^^-i:u.ir)dri^rj^^)uMdy' 

^" h{\\viA\)v, + b',{\\nA\){!^v{y)dy + v,)u.,]) 

+ (feo(||u,||)fe', (^^^^) {^^^^^ - ^^^^ a'".(y)d, + .o)) -.o) ' 
which is a bounded linear operator on X . 
Theorem 4.7. // 

y.^ inf {m(s)}> ||F^J|, (4.25) 

.•5G[0.m] " 

then the non-trivial steady state u» is locally asymptotically stable. 

Proof. Note that A — M has compact resolvent (where M is the mortality operator 
introduced in p.l9p ). see the proof of Lemma lOl and generates a positive and irreducible 
semigroup, see the proof of Lemma [3.21 Therefore its point spectrum ap{A — M) is not 
empty, see [H C-III Theorem 3.7]. We also noted before that it generates a contraction 
semigroup. In fact, the semigroup 7i(t) generated hy A — M satisfies ||7i(t)||;^. = 1 
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for every t > 0, hence = ujo{A — M) = s{A — M). On the other hand we have 
I |exp {tF^^ } 1 1 < exp I I^F^^ 1 1 }, hence for the growth bound of the semigroup T2{t) 
generated by we have wq (^'u.) — | |-^u. 1 1- We also note that M + z/J is a dissipative 
operator and it generates a contraction semigroup, hence for the growth bound of the 
semigroup Tz{t) generated by M we have ujo^M) < —u, and ||73(i)|| < exp{— i^t}. 
FinaUy, by applying a version of the Trotter product formula (see e.g. |13l Corollary 
5.8]) we obtain a;o(A - M + F^^ + M) = iuo(A + F^ J < which shows that the non- 
trivial steady state is locally asymptotically stable. □ 

Remark 4.8 The stability condition in Theorem 14.71 mav seem very restrictive at the 
first glance. It may however be the case that the norm of the linearisation jlF^^H is 
small, especially since for our model we naturally have /3g, 6q, /32, &2 — ^- On the other 
hand, even natural mortality (note that is natural mortality combined with infection 
induced mortality) is very high for example in the case of the mosquito Aedes aegypti, 
which is one of the target species for the introduction of Wolbachia, see [31]. We also 
note that the operator —F^^ may well be positive if the conditions of Theorem 13.51 are 
satisfied, but there are no readily available results concerning the question that how does 
the infimum of the spectrum of an operator changes under positive perturbations. On 
the other hand if F^^ is dissipative and ly > then similar arguments as used in the 
proof of Theorem 14.71 show that the steady state is locally asymptotically stable. 

Remark 4.9 As we pointed out earlier, it is not possible to talk about linear stabi- 
lity of the trivial steady state in the usual sense, as the nonlinearity F is not Frechet 
differentiable at 0. We established however global existence of solutions for any initial 
condition, and since condition (ii) in Theorem [33] requires that the spectral bound of the 
linear part is positive in a small neighbourhood of it may be intuitively plausible to 
expect (for example utilising the variation formula (|2.16p ) that the trivial steady state 
is unstable and the positive steady state is actually globally asymptotically stable if the 
condition of Theorem 14.71 holds true. 

5. Concluding remarks 

In this paper we introduced and analysed a nonlinear structured population model with 
diffusion in the state space. Individuals in the population are structured with respect 
to infection (for example bacterium) load, hence we used Wentzell boundary condition 
at the uninfected state x = 0. Our model primarily intended to describe the evolution 
of Wolbachia infection in an arthropod, for example mosquito population. Wolbachia 
is a reproductive parasite and it affects the reproductive mechanisms of its host in an 
intriguing fashion. Here we focused on a cytoplasmic incompatibility (CI) inducing strain. 
Following [6] [25] we adopted the view of partial CI, namely, that a female can produce 
viable offspring only when mating with a male who has lower infection load. Therefore 
the functions /?2,&2 are assumed to be monotone decreasing. Note however, that these 
assumptions are not necessary to establish any of our results, such as Theorem 12.11 and 
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Theorem 13.51 In fact, in our main result Theorem 13.51 the crucial assumption is the 
one which concerns the strict monotonocity of the functions /3o;^o- The necessity of 
condition (i) in Theorem 13.51 is also in agreement with the fact that CI itself does not 
regulate population growth, it only provides infected individuals with a reproductive 
advantage, see e.g. [T3J[TC1[31]. The main difficulty in the mathematical analysis of our 
model also arises from this crucial assumption of the density dependent CI. In particular, 
physiologically structured population models with infinite dimensional nonlinearities and 
with distributed recruitment terms, such as the one used here, are usually notoriously 
difficult to analyse, see e.g. [SI [H]- The special form of the nonlinearity naturally 
appearing in P2 and 62 also implies that although the nonlinear operator F is locally 
Lipschitz continuous at it is not Frechet differentiable. This also implies that even 
though intuitively one may expect that if the conditions of Theorem 13.51 and Theorem 
14.71 hold true then the trivial steady state is unstable, we cannot discuss linear stability 
of the trivial steady state using standard stability results from scmilinear theory as 
developed for example in [2E1 EE] ■ 
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